Nuclear Spin Effects in Optical Lattice Clocks 



o 
o 

(N 



Martin M. Boyd, Tanya Zelevinsky, Andrew D. Ludlow, Sebastian 
Blatt, Thomas Zanon-Willette, Seth M. Foreman, and Jun Ye 
JILA, National Institute of Standards and Technology and University of Colorado, 
Department of Physics, University of Colorado, Boulder, CO 80309-0440 
(Dated: February 1, 2008) 

We present a detailed experimental and theoretical study of the effect of nuclear spin on the 
performance of optical lattice clocks. With a state-mixing theory including spin-orbit and hyperfine 
interactions, we describe the origin of the 1 So- 3 Po clock transition and the differential g-factor be- 
tween the two clock states for alkaline-earth (-like) atoms, using 87 Sr as an example. Clock frequency 
shifts due to magnetic and optical fields are discussed with an emphasis on those relating to nuclear 
structure. An experimental determination of the differential g-factor in 87 Sr is performed and is 
in good agreement with theory. The magnitude of the tensor light shift on the clock states is also 
explored experimentally. State specific measurements with controlled nuclear spin polarization are 
discussed as a method to reduce the nuclear spin-related systematic effects to below 10 -17 in lattice 
clocks. 
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Optical clocks [lj based on alkaline-earth atoms con- 
fined in an optical lattice [1J are being intensively ex- 
plored as a route to improve state of the art clock accu- 
racy and precision. Pursuit of such clocks is motivated 
mainly by the benefits of Lamb-Dicke confinement which 
allows high spectral resolution [H, 0] , and high accuracy 
[H, fl 0, H| with the suppression of motional effects, while 
the impact of the lattice potential can be eliminated using 
the Stark cancelation technique 0, [T3, EH, E3] • Lattice 
clocks have the potential to reach the impressive accu- 
racy level of trapped ion systems, such as the Hg + opti- 
cal clock jl3(, while having an improved stability due to 
the large number of atoms involved in the measurement. 
Most of the work performed thus far for lattice clocks has 
been focused on the nuclear-spin induced 1 Sq- 3 Pq tran- 
sition in 87 Sr. Recent experimental results are promis- 
ing for development of lattice clocks as high performance 
optical frequency standards. These include the confir- 
mation that hyperpolarizability effects will not limit the 
clock accuracy at the 10 -17 level [l^, observation of tran- 
sition resonances as narrow as 1.8 Hz Q, and the excel- 
lent agreement between high accuracy frequency mea- 
surements performed by three independent laboratories 
[E IE 0, Hi with clock systematics associated with the lat- 
tice technique now controlled below 10 -15 @. A main 
effort of the recent accuracy evaluations has been to min- 
imize the effect that nuclear spin (1 = 9/2 for 87 Sr) has 
on the performance of the clock. Specifically, a linear 
Zeeman shift is present due to the same hyperfine inter- 
action which provides the clock transition, and magnetic 
sublevel-dependent light shifts exist, which can compli- 
cate the stark cancelation techniques. To reach accuracy 
levels below 10 -17 , these effects need to be characterized 
and controlled. 

The long coherence time of the clock states in alkaline 
earth atoms also makes the lattice clock an intriguing 
system for quantum information processing. The closed 
electronic shell should allow independent control of elec- 
tronic and nuclear angular momenta, as well as protec- 
tion of the nuclear spin from environmental perturbation, 



providing a robust system for coherent manipulation [Tj]. 
Recently, protocols have been presented for entangling 
nuclear spins in these systems using cold collisions [151 ] 
and performing coherent nuclear spin operations while 
cooling the system via the electronic transition [l6| . 

Precise characterization of the effects of electronic and 
nuclear angular-momentum-interactions and the resul- 
tant state mixing is essential to lattice clocks and po- 
tential quantum information experiments, and therefore 
is the central focus of this work. The organization of 
this paper is as follows. First, state mixing is discussed 
in terms of the origin of the clock transition as well as 
a basis for evaluating external field sensitivities on the 
clock transition. In the next two sections, nuclear-spin 
related shifts of the clock states due to both magnetic 
fields and the lattice trapping potential are discussed. 
The theoretical development is presented for a general 
alkaline-earth type structure, using 87 Sr only as an ex- 
ample (Fig. [T]) , so that the results can be applied to other 
species with similar level structure, such as Mg, Ca, Yb, 
Hg, Zn, Cd, Al + , and In + . Following the theoretical dis- 
cussion is a detailed experimental investigation of these 
nuclear spin related effects in 87 Sr, and a comparison to 
the theory sections. Finally, the results are discussed in 
the context of the performance of optical lattice clocks, 
including a comparison with recent proposals to induce 
the clock transition using external fields in order to elim- 
inate nuclear spin effects [13, H GJ, HE SH, S2|. The 
appendix contains additional details on the state mixing 
and magnetic sensitivity calculations. 



STATE MIXING IN THE nsnp 
CONFIGURATION 



To describe the two-electron system in intermediate 
coupling, we follow the method of Breit and Wills (23j 
and Lurio [24| and write the four real states of the ns np 
configuration as expansions of pure spin-orbit (LS) cou- 
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FIG. 1: (color online) Simplified 87 Sr energy level diagram 
(not to scale). Relevant optical transitions discussed in the 
text are shown as solid arrows, with corresponding wave- 
lengths given in nanometers. Hyperfine structure sublevels 
are labeled by total angular momentum F, and the magnetic 
dipole (A) and electric quadrupole (Q, equivalent to the hy- 
perfine B coefficient) coupling constants are listed in the inset. 
State mixing of the x Pi and 3 Pi states due to the spin-orbit 
interaction is shown as a dashed arrow. Dotted arrows repre- 
sent the hyperfine induced state mixing of the 3 Po state with 
the other F — 9/2 states in the 5s5p manifold. 



states having the same total spin F, mixing the pure 3 Po 
state with the 3 Pi, 3 P2 and x Pi states. 

I 3 i%> = | 3 P °> +ao| 3 Pi) +/3 rPi) +7o| 3 P 2 °). (2) 

The HFI mixing coefficients ao, 0o, and 70 (2x 10~ 4 , — 4x 
1CT 6 , and 4 x 10" 6 respectively for 87 Sr) are defined in 
Ea.[T51of the appendix and can be related to the hyperfine 
splitting in the P states, the fine structure splitting in the 
3 P states, and the coupling coefficients a and j3 [23. |2^|. 
The 3 Po state can also be written as a combination of 
pure states using Eq. [T] 



| 3 Po)=| 3 ^o> + («o«-/3o/3)| 3 Pi°) 

+ (Qo/3 + /3oa)| 1 Pi>+7o| 3 P 2 ). 



(3) 



The HFI mixing enables a non-zero electric-dipole tran- 
sition via the pure 1 P® state, with a lifetime which can 
be calculated given the spin-orbit and HFI mixing coef- 
ficients, the 3 Pi lifetime, and the wavelengths (A) of the 
3 Po and 3 Pi transitions from the ground state [391 ] . 
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Here the intermediate coupling coefficients a and f3 
(0.9996 and -0.0286 respectively for Sr) represent the 
strength of the spin-orbit induced state mixing between 
singlet and triplet levels, and can be determined from 
experimentally measured lifetimes of *Pi and 3 Pi (see 
Eq. [15] in the appendix) . This mixing process results in 
a weakly allowed 1 Sq- 5 P\ transition (which would other- 
wise be spin- forbidden), and has been used for a variety 
of experiments spanning different fields of atomic physics. 
In recent years, these intercombination transitions have 
provided a unique testing ground for studies of narrow- 
line cooling in Sr [H, Sa EE HI HI and Ca [HHH, as 
well as the previously unexplored regime of photoassocia- 
tion using long lived states [H, [H, [3J] . These transitions 
have also received considerable attention as potential op- 
tical frequency standards [35|, HH, [37j], owing mainly to 
the high line quality factors and insensitivity to external 
fields. Fundamental symmetry measurements, relevant 
to searches of physics beyond the standard model, have 
also made use of this transition in Hg [38J. Furthermore, 
the lack of hyperfine structure in the bosonic isotopes 
(I = 0) can simplify comparison between experiment and 
theory. 

The hyperfine interaction (HFI) in fermionic isotopes 
provides an additional state mixing mechanism between 



In the case of Sr, the result is a natural lifetime on the 
order of 100 seconds HO, El|, compared to that of a 
bosonic isotope where the lifetime approaches 1000 years 
[4ll ] . Although the 100 second coherence time of the 
excited state exceeds other practical limitations in cur- 
rent experiments, such as laser stability or lattice life- 
time, coherence times approaching one second have been 
achieved Q. The high spectral resolution has allowed a 
study of nuclear-spin related effects in the lattice clock 
system discussed below. 

The level structure and state mixing discussed here 
are summarized in a simplified energy diagram, shown 
in Fig. [TJ which gives the relevant atomic structure and 
optical transitions for the 5s5p configuration in Sr. 



II. THE EFFECT OF EXTERNAL MAGNETIC 
FIELDS 

With the obvious advantages in spectroscopic precision 
of the 1 S'o- 3 Po transition in an optical lattice, the sensi- 
tivity of the clock transition to external field shifts is a 
central issue in developing the lattice clock as an atomic 
frequency standard. To evaluate the magnetic sensitivity 
of the clock states, we follow the treatment of Ref. [24| for 
the intermediate coupling regime described by Eqns.QJH] 
in the presence of a weak magnetic field. A more general 
treatment for the case of intermediate fields is provided 
in the appendix. The Hamiltonian for the Zccman inter- 
action in the presence of a weak magnetic field B along 
the z-axis is given as 



H z = {g s S z + giL z - g!l z )fx B. 



(5) 



3 



Here g s ~ 2 and gi = 1 are the spin and orbital an- 
gular momentum <?-factors, and S^, L Z1 and J a are the 
z-components of the electron spin, orbital, and nuclear 
spin angular momentum respectively. The nuclear g- 
factor, gi, is given by gi= ^~}t^ > where /i/ is the nuclear 
magnetic moment, ad is the diamagnetic correction and 
/io =i 7f-- Here, \ib is the Bohr magneton, and h is Planck's 
constant. For 87 Sr, the nuclear magnetic momement and 
diamagnetic correction are \ij = — 1.0924(7)/zat [H| and 
cr d = 0.00345 [43[ respectively, where fi^ is the nuclear 
magneton. In the absence of state mixing, the 3 Pq g- 
factor would be identical to the 1 Sq ^-factor (assuming 
the diamagnetic effect differs by a negligible amount for 
different electronic states), equal to gi. However since 
the HFI modifies the 3 Pq wavcfunction, a differential g- 
factor, dg, exists between the two states. This can be 
interpreted as a paramagnetic shift arising due to the 
distortion of the electronic orbitals in the triplet state, 
and hence the magnetic moment [44 |. Sg is given by 
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( 3 P \H Z \ 3 P ) - { s P8\H z \ 3 P°) 
rriFfAoB 

- , „ m ( 3 Plm F \H z \ 3 P[\F = I,m F ) 
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Using the matrix element given in the appendix for 
87 Sr (I = 9/2), we find { 3 P^m F \H z \ 3 P^F = l,m F )= 

rriFHoB, corresponding to a modification of the 3 Pq 
g-factor by ~60%. Note that the sign in Eq. [6] differs 
from that reported in [39|, |44| due to our choice of sign 
for the nuclear term in the Zeeman Hamiltonian (oppo- 
site of that found in Ref. 24}). The resulting linear 
Zeeman shift — —SgmF/MiB of the 1 Sq- 3 Pq transition 
is on the order of ~110xm F Hz/G (1 G = 10~ 4 Tesla). 
This is an important effect for the development of lattice 
clocks, as stray magnetic fields can broaden the clock 
transition (deteriorate the stability) if multiple sublevels 
are used. Furthermore, imbalanced population among 
the sublevels or mixed probe polarizations can cause fre- 
quency errors due to line shape asymmetries or shifts. 
It has been demonstrated that if a narrow resonance is 
achieved (10 Hz in the case of Ref. 0), these systematics 
can be controlled at 5xl0 -16 for stray fields of less than 
5 mG. To reduce this effect, one could employ narrower 
resonances or magnetic shielding. 

An alternative measurement scheme is to measure 
the average transition frequency between itif and —wif 
states of to cancel the frequency shifts. This requires 
application of a bias field to resolve the sublevels, and 
therefore the second order Zeeman shift A^' must be 
considered. The two clock states are both J = so the 
shift Ag arises from levels separated in energy by the 
fine-structure splitting, as opposed to the more tradi- 
tional case of alkali (-like) atoms where the second order 
shift arises from nearby hyperfine levels. The shift of 
the clock transition is dominated by the interaction of 
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FIG. 2: (color online) A Breit-Rabi diagram for the 1 So- 3 Po 
clock transition using Eq. [22] with Sg^io = —109 Hz/G. Inset 
shows the linear nature of the clock shifts at the fields relevant 
for the measurement described in the text. 

the 3 Po and 3 Pi states since the ground state is sepa- 
rated from all other energy levels by optical frequencies. 
Therefore, the total Zeeman shift of the clock transition 
As is given by 
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The frequency difference in the denominator is mainly 
due to the fine-structure splitting and is nearly indepen- 
dent of F', and can therefore be pulled out of the sum- 
mation. In terms of the pure states, and ignoring terms 
of order ao, /3q, /3 2 , and smaller, we have 



,(2) 
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where we have used the matrix elements given in the 
appendix for the case F = 9/2. From Eq. [5] the sec- 
ond order Zeeman shift (given in Hz for a magnetic field 
given in Gauss) for 87 Sr is A^ 2) =-0.233S 2 . This is con- 
sistent with the results obtained in Ref. (2(| and [45[ for 
the bosonic isotope. Inclusion of the hyperfine splitting 
into the frequency difference in the denominator of Eq. [7| 
yields an additional term in the second order shift pro- 
portional to m F which is more that 10 -6 times smaller 
than the main effect, and therefore negligible. Notably, 
the fractional frequency shift due to the second order 
Zeeman effect of 5xl0 -16 G~ 2 is nearly 10 8 times smaller 
than that of the Cs 0, |4?| clock transition, and more 
than an order of mag nitude smaller than that present in 
Hg+ [HI, Sr+ m,|lf,andYb + |o|,[5l| ion optical clocks. 

A Breit-Rabi like diagram is shown in Fig. giving 
the shift of the ^^Sq^Pq transition frequency for different 
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nip sublevels (assuming Am = for n transitions), as a 
function of magnetic field. The calculation is performed 
using an analytical Breit-Rabi formula (Eq. |22|) provided 
in the appendix. The result is indistinguishable from the 
perturbative derivation in this section, even for fields as 
large as 10 4 G. 

III. THE EFFECT OF THE OPTICAL LATTICE 
POTENTIAL 

In this section we consider the effect of the confining 
potential on the energy shifts of the nuclear sublevels. In 
the presence of a lattice potential of depth Up, formed by 
a laser linearly polarized along the axis of quantization 
defined by an external magnetic field B, the level shift of 
a clock state (hA g / e ) from its bare energy is given by 

A e = - m F (g I + Sg)^i B - nfj^- - K^£m F ^- 

Er tjR 

-kZ (3m 2 F -F(F + l))^- 

U U W 

A n s Ut V,. Ut 

Ag = - m Fgi ^ B -Kg— Kg ^ F — 

Er Er 

-K T a (3m|-F(F + l))g. 

Here, k s , k v , and k t are shift coefficients proportional 
to the scalar, vector (or axial), and tensor polarizabil- 
ities, and subscripts e and g refer to the excited ( 3 Po) 
and ground ^Sq) states respectively. Er is the energy 
of a lattice photon recoil and Ut/Er characterizes the 
lattice intensity. The vector (cx mp) and tensor (cx mp) 
light shift terms arise solely from the nuclear structure 
and depend on the orientation of the light polarization 
and the bias magnetic field. The tensor shift coefficient 
includes a geometric scaling factor which varies with the 
relative angle </> of the laser polarization axis and the axis 
of quantization, as 3cos 2 <fi — 1. The vector shift, which 
can be described as an pseudo-magnetic field along the 
propagation axis of the trapping laser, depends on the 
trapping geometry in two ways. First, the size of the 
effect is scaled by the degree of elliptical polarization £, 
where £ = (£ = ±1) represents perfect linear (circular) 
polarization. Second, for the situation described here, 
the effect of the vector light shift is expected to be orders 
of magnitude smaller than the Zeeman effect, justifying 
the use of the bias magnetic field direction as the quan- 
tization axis for all of the mp terms in Eq. [TJJ Hence 
the shift coefficient depends on the relative angle be- 
tween the pseudo-magnetic and the bias magnetic fields, 
vanishing in the case of orthogonal orientation A 
more general description of the tensor and vector effects 
in alkaline-earth systems for the case of arbitrary ellipti- 
cal polarization can be found in Ref. [lj| • Calculations of 
the scalar, vector, and tensor shift coefficients have been 
performed elsewhere for Sr, Yb, and Hg Emm 
and will not be discussed here. Hyperpolarizability ef- 
fects (cx U T ) [I, [M ED, O are ignored in Eq. M as they 



are negligible in 87 Sr at the level of 10~ 17 for the range of 
lattice intensities used in current experiments [T3 |. The 
second order Zeeman term has been omitted but is also 
present. 

Using Eq. [9] we can write the frequency of a n- 
transition (Amp = 0) from a ground state mj? as 

f* ra , = ~ (Ak S - Ak t F(F + 1)) ^ 

- (WWc + A^Sm 2 ^^ ( 10 ) 

V / Er 

— 8gm F ^j,oB, 

where the shift coefficients due to the differential polar- 
izabilities are represented as Ak, and v c is the bare clock 
frequency. The basic principle of the lattice clock tech- 
nique is to tune the lattice wavelength (and hence the 
polarizabilities) such that the intensity-dependent fre- 
quency shift terms are reduced to zero. Due to the mp- 
dependence of the third term of Eq. [TUJ the Stark shifts 
cannot be completely compensated for all of the sublevels 
simultaneously. Or equivalently, the magic wavelength 
will be different depending on the sublevel used. The 
significance of this effect depends on the magnitude of 
the tensor and vector terms. Fortunately, in the case of 
the 1 5'o- 3 Po transition the clock states are nearly scalar, 
and hence these effects are expected to be quite small. 
While theoretical estimates for the polarizabilities have 
been made, experimental measurements are unavailable 
for the vector and tensor terms. The frequencies of o 
(Amp = ±1) transitions from a ground state are 
similar to the 7r-transitions, given by 

v ± = v c - (Ak s - Ak t F(F +1])^- 

CTm F V / Er, 

-(iKY(m F ±l)-Kjm F )^j^ 

— (Ke3(m F ± l) 2 — K^Sm'p) 

V / Er 

- (±gi + 8g(m F ± l))fi, B. 

IV. EXPERIMENTAL DETERMINATION OF 
FIELD SENSITIVITIES 

To explore the magnitude of the various rnp-dependent 
shifts in Eq. 1101 a differential measurement scheme can 
be used to eliminate the large shifts common to all levels. 
Using resolved sublevels one can extract m_p sensitivities 
by measuring the splitting of neighboring states. This is 
the approach taken here. A diagram of our spectroscopic 
setup is shown in Fig. [3ja). 87 Sr atoms are captured 
from a thermal beam into a magneto-optical trap (MOT), 
based on the 1 Sq- P\ cycling transition. The atoms are 
then transferred to a second stage MOT for narrow line 
cooling using a dual frequency technique (2q |. Full de- 
tails of the cooling and trapping system used in this 
work are discussed elsewhere [a, [28[ . During the cooling 
process, a vertical one-dimensional lattice is overlapped 
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FIG. 3: (color online) (a) Schematic of the experimental ap- 
paratus used here. Atoms are confined in a nearly vertical 
optical lattice formed by a retro-reflected 813 nm laser. A 
698 nm probe laser is co-aligned with the lattice. The probe 
polarization Ep can be varied by an angle 6 relative to that of 
the linear lattice polarization El- A pair of Helmholtz coils 
(blue) is used to apply a magnetic field along the lattice po- 
larization axis, (b) Nuclear structure of the 1 So and 3 Po clock 
states. The large nuclear spin (J = 9/2) results in 28 total 
transitions, and the labels ir, a + , and a~ represent transi- 
tions where rar changes by 0, +1, and —1 respectively, (c) 
Observation of the clock transition without a bias magnetic 
field. The 3 Po population (in arbitrary units) is plotted (blue 
dots) versus the probe laser frequency for 9 — 0, and a fit to 
a sinc-squared lineshape yields a Fourier-limited linewidth of 
10.7(3) Hz. Linewidths as narrow as 5 Hz have been observed 
under similar conditions and when the probe time is extended 
to 500 ms. 



with the atom cloud. We typically load — 10 4 atoms into 
the lattice at a temperature of — 1.5/iK. The lattice is 
operated at the Stark cancelation wavelength 0, [13] of 
813.4280(5) nm with a trap depth of U = 35E R . A 
Helmholtz coil pair provides a field along the lattice po- 
larization axis for resolved sub-level spectroscopy. Two 
other coil pairs are used along the other axes to zero the 
orthogonal fields. The spectroscopy sequence for the 
3 Po clock transition begins with an 80 ms Rabi pulse from 
a highly stabilized diode laser [53| that is co-propagated 
with the lattice laser. The polarization of the probe laser 
is linear at an angle 8 relative to that of the lattice. A 
shelved detection scheme is used, where the ground state 
population is measured using the 1 Sq- 1 Pi transition. The 
3 Po population is then measured by pumping the atoms 
through intermediate states using 3 Po- 3 Si, 3 P2- 3 'Si, and 
the natural decay of 3 Pi , before applying a second 1 S'o- 
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FIG. 4: (color online) Observation of the 1 S l - 3 Po ir- 
transitions (6 = 0) in the presence of a 0.58 G magnetic field. 
Data is shown in grey and a fit to the eight observable line- 
shapes is shown as a blue curve. The peaks are labeled by 
the ground state mF-sublevel of the transition. The relative 
transition amplitudes for the different sublevels are strongly 
influenced by the Clebsch-Gordan coefficients. Here, transi- 
tion linewidths of 10 Hz are used. Spectra as narrow as 1.8 
Hz can be achieved under similar conditions if the probe time 
is extended to 500 ms. 



1 Pi pulse. The 461 nm pulse is destructive, so for each 
frequency step of the probe laser the ~800 ms loading 
and cooling cycle is repeated. 

When 7r polarization is used for spectroscopy (8 = 0), 
the large nuclear spin provides ten possible transitions, 
as shown schematically in Fig. [3£b). Figure [3Jc) shows 
a spectroscopic measurement of these states in the ab- 
sence of a bias magnetic field. The suppression of mo- 
tional effects provided by the lattice confinement allows 
observation of extremely narrow lines [H, 0, [n| , in this 
case having Fourier-limited full width at half maximum 
(FWHM) of -10 Hz (quality factor of 4 x 10 13 ). In our 
current apparatus the linewidth limitation is 5 Hz with 
degenerate sublevels and 1.8 Hz when the degeneracy is 
removed [3j. The high spectral resolution allows for the 
study of nuclear spin effects at small bias fields, as the 
ten sublevels can easily be resolved with a few hundred 
mG. An example of this is shown in Fig.0] where the ten 
transitions are observed in the presence of a 0.58 G bias 
field. This is important for achieving a high accuracy 
measurement of Sg as the contribution from magnetic- 
field-induced state mixing is negligible. To extract the 
desired shift coefficients we note that for the tt transi- 
tions we have a frequency gap between neighboring lines 
of 



v U T t Ut (12) 

-5gu B - Ak v £ — - A/t T 3 2m F - 1) — . 
yP ^E R y ' Er 



From Eq. 1121 we see that by measuring the differences in 
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FIG. 5: (color online) Observation of the 18 a transitions 
when the probe laser polarization is orthogonal to that of the 
lattice (9 = -|). Here, a field of 0.69 G is used. The spectro- 
scopic data is shown in grey and a fit to the data is shown 
as a blue curve. Peak labels give the ground state sublevel of 
the transition, as well as the excitation polarization. 



frequency of two spectroscopic features, the three terms 
of interest (Sg, Ak v , and Ak t ) can be determined inde- 
pendently. The differential g factor can be determined 
by varying the magnetic field. The contribution of the 
last two terms can be extracted by varying the inten- 
sity of the standing wave trap, and can be independently 
determined since only the tensor shift depends on mp. 

While the ir transitions allow a simple determination 
of Sg, the measurement requires a careful calibration of 
the magnetic field and a precise control of the probe 
laser frequency over the ~500 seconds required to pro- 
duce a scan such as in Fig. |U Any linear laser drift 
will appear in the form of a smaller or larger Sg, de- 
pending on the laser scan direction. Furthermore, the 
measurement can not be used to determine the sign of 
Sg as an opposite sign would yield an identical spectral 
pattern. In an alternative measurement scheme, we in- 
stead polarize the probe laser perpendicular to the lattice 
polarization (0 — -|) to excite both er + and o~ tran- 
sitions. In this configuration, 18 spectral features are 
observed and easily identified (Fig. [5]). Ignoring small 
shifts due to the lattice potential, Sg is given by extract- 
ing the frequency splitting between adjacent transitions 
of a given polarization (all cr + or all a~ transitions) as 
/<r±.m F — u a ± ~ v a ± =—Sgfi B. If we also measure the 

frequency difference between a + and o~ transitions from 
the same sublevel, fd,m F =v + — v - =—2(gi + 8g)goB, 

we find that the differential g-factor can be determined 
from the ratio of these frequencies as 
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FIG. 6: (color online) Calculation of the 18 a transition fre- 
quencies in the presence of a 1 G bias field, including the influ- 
ence of Clebsch-Gordan coefficients. The green (red) curves 
show the a + (cr _ ) transitions, (a) Spectral pattern for g- 
factors giUo = —185 Hz/G and Sg/io — —109 Hz/G. (b) Same 
pattern as in (a) but with Sggo = +109 Hz/G. The qualita- 
tive difference in the relative positions of the transitions allows 
determination of the sign of Sg compared to that of gi. 



In this case, prior knowledge of the magnetic field is not 
required for the evaluation, nor is a series of measure- 
ment at different fields, as Sg is instead directly deter- 
mined from the line splitting and the known 1 Sq g factor 
gi. The field calibration and the Sg measurement are in 
fact done simultaneously, making the method immune to 
some systematics which could mimic a false field, such as 
linear laser drift during a spectroscopic scan or slow mag- 
netic field variations. Using the a transitions also elim- 
inates the sign ambiguity which persists when using the 
7r transitions for measuring Sg. While we can not extract 
the absolute sign, the recovered spectrum is sensitive to 
the relative sign between gi and Sg. This is shown explic- 
itly in Fig. [5] where the positions of the transitions have 
been calculated in the presence of a ~1 G magnetic field. 
Figure [6]^a) shows the spectrum when the signs of gi and 
Sg are the same while in Fig. [HJb) the signs are oppo- 
site. The two plots show a qualitative difference between 
the two possible cases. Comparing Fig. [5] and Fig. [6] it is 
obvious that the hyperfine interaction increases the mag- 
nitude of the 3 Po (/-factor (Sg has the same sign as gi). 
We state this point explicitly because of recent inconsis- 
tencies in theoretical estimates of the relative sign of Sg 
and gi in the 87 Sr literature 

To extract the magnitude of Sg, data such as in Fig. O 
are fit with eighteen Lorentzian lines, and the relevant 
splitting frequencies fd,m F and f a ± are extracted. Due 
to the large number of spectral features, each experimen- 
tal spectrum yields 16 measurements of Sg. A total of 
31 full spectra was taken, resulting in an average value 
of Sggo = —108.4(4) Hz/G where the uncertainty is the 
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FIG. 7: (color online) Summary of ^-measurements for dif- 
ferent lattice intensities. Each data point (and uncertainty) 
represents the Sg value extracted from a full a spectrum 
such as in Fig. [5] Linear extrapolation (red line) to zero lat- 
tice intensity yields a value —108.4(1) Hz/G. 



standard deviation of the measured value. To check for 
sources of systematic error, the magnetic field was varied 
to confirm the field independence of the measurement. 
We also varied the clock laser intensity by an order of 
magnitude to check for Stark and line pulling effects. It 
is also necessary to consider potential measurement er- 
rors due to the optical lattice since in general the splitting 
frequencies fd,m F and f a ± will depend on the vector and 
tensor light shifts. For fixed fields, the vector shift is in- 
distinguishable from the linear Zeeman shift (see Eas.llOt- 
[T2")) and can lead to errors in calibrating the field for a Sg 
measurement. In this work, a high quality linear polar- 
izer (10 -4 ) is used which would in principle eliminate the 
vector shift. The nearly orthogonal orientation should 
further reduce the shift. However, any birefringence of 
the vacuum windows or misalignment between the lattice 
polarization axis and the magnetic field axis can lead to 
a non-zero value of the vector shift. To measure this ef- 
fect in our system, we varied the trapping depth over a 
range of ~ (0.6 - 1.7)f7o and extrapolated Sg to zero in- 
tensity, as shown in Fig. [71 Note that this measurement 
also checks for possible errors due to scalar and tensor 
polarizabilites as their effects also scale linearly with the 
trap intensity. We found that the (^-measurement was 
affected by the lattice potential by less then 0.1%, well 
below the uncertainty quoted above. 

Unlike the vector shift, the tensor contribution to the 
sublevel splitting is distinguishable from the magnetic 
contribution even for fixed fields. Adjacent a transitions 
can be used to measure An T and due to the m 2 F de- 
pendence of the tensor shift. An appropriate choice of 
transition comparisons results in a measurement of the 
tensor shift without any contributions from magnetic or 
vector terms. To enhance the sensitivity of our measure- 
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FIG. 8: (color online) Measurement of the tensor shift coef- 
ficients Ak t (blue triangles), and (green circles), using a 
spectra and Eq. 1141 The measured coefficients show no sta- 
tistically significant trap depth dependence while varying the 
depth from 0.85-1.7 U . 



ment we focus mainly on the transitions originating from 
states with large rriF', for example, we find that 



Ak 1 



/, 



<T+,m J 7 = 7/2 



-7/2 



42^21 



=7/2 



-7/2 



(14) 



84^ 



while similar combinations can be used to isolate the dif- 
ferential tensor shift from the a~ data as well as the 
tensor shift coefficient of the 1 So state. From the a split- 
ting data we find Ak t = 0.03(8) Hz/f/ and |k^|=0.02(4) 
Hz/[/o- The data for these measurements is shown in 
Fig. [H Similarly, we extracted the tensor shift coeffi- 
cient from 7r spectra, exploiting the m_F-dependent term 
in Eq. [TJ yielding Ak t = 0.02(7) Hz/(7 - The measure- 
ments here are consistent with zero and were not found 
to depend on the trapping depth used for a range of 0.85- 
1.7 £/o, and hence are interpreted as conservative upper 
limits to the shift coefficients. The error bars represent 
the standard deviation of many measurements, with the 
scatter in the data due mainly to laser frequency noise 
and slight under sampling of the peaks. It is worth noting 
that the tensor shift of the clock transition is expected to 
be dominated by the 3 Po shift, and therefore, the limit 
on can be used as an additional estimate for the up- 
per limit on Ak t . Improvements on these limits can be 
made by going to larger trap intensities to enhance sen- 
sitivity, as well as by directly stabilizing the clock laser 
to components of interest for improved averaging. 

Table [T] summarizes the measured sensitivities to mag- 
netic fields and the lattice potential. The Stark shift 
coefficients for linear polarization at 813.4280(5) nm are 
given in units of Hz /{Ut I Er). For completeness, a recent 
measurement of the second order Zeeman shift using 88 Sr 
has been included [45[ , as well as the measured shift coef- 
ficient A7 for the hyperpolarizability [12j and the upper 
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TABLE I: Measured Field Sensitivities for 87 Sr 



Sensitivity 


Value 


Units 


Ref. 


Al!' IravB 
Af/B 2 
Ak t 
Ak t 

K e 
K 

A-y 


-108.4(4) 
-0.233(5) 
6(20) xlO" 4 
9(23) xlO" 4 
5(10) xlO" 4 
-3(7) xlO" 3 
7(6)xl0~ 6 


Hz/G 

Hz/G 2 
Hz/(U T /E R ) 
Uz/(Ut/E r ) 
Hz/(U T /E R ) 
Hz/(U t /Er) 
Rz/(Ut/E r ) 2 


This work 

[45j a 
This work 6 
This work c 
This work c 

[6] d 

[12] d 


a Measured for 88 Sr 

b Measured with n spectra 

c Measured with a ± spectra 

d Measured with degenerate sublevels 



limit for the overall linear lattice shift coefficient k from 
our recent clock measurement [6j. While we were able 
to confirm that the vector shift effect is small and con- 
sistent with zero in our system, we do not report a limit 
for the vector shift coefficient Ak v due to uncertainty 
in the lattice polarization purity and orientation relative 
to the quantization axis. In future measurements, use of 
circular trap polarization can enhance the measurement 
precision of An v by at least two orders of magnitude. 

Although only upper limits are reported here, the re- 
sult can be used to estimate accuracy and linewidth lim- 
itations for lattice clocks. For example, in the absence 
of magnetic fields, the tensor shift can cause line broad- 
ening of the transition for unpolarized samples. Given 
the transition amplitudes in Fig. 01 the upper limit for 
line broadening, derived from the tensor shift coefficients 
discussed above, is 5 Hz at Uo- The tensor shift also 
results in a different magic wavelength for different mj? 
sublevels, which is constrained here to the few picometer 
level. 



V. COMPARISON OF THE Sg MEASUREMENT 
WITH THEORY AND 3 P LIFETIME ESTIMATE 

The precise measurement of Sg provides an opportu- 
nity to compare various atomic hyperfine interaction the- 
ories to the experiment. To calculate the mixing param- 
eters ao and f3o (defined in Eq. [T|5] of the Appendix), 
we first try the simplest approach using the standard 
Breit- Wills (BW) theory [H, HI to relate the mixing 
parameters to the measured triplet hyperfine splitting 
(hfs). The parameters a (0.9996) and (5 (-0.0286(3)) 
are calculated from recent determinations of the 3 Pi [32j 
and 1 P\ [54J lifetimes. The relevant singlet and triplet 
sing le-electron hyperfine coefficients are taken from Ref. 
[5a |. From this calculation we find ao = 2.37(1) x 10 -4 , 
A) = -4.12(1) x 10~ 6 , and 70 = 4.72(1) x 10~ 6 , resulting 
in Sgfio = —109.1(1) Hz/G . Using the mixing values in 
conjunction with Eq. [4] we find that the 3 Pq lifetime is 
152(2) s. The agreement with the measured g-factor is 
excellent, however the BW-theory is known to have prob- 
lems predicting the 1 P\ characteristics based on those of 



the triplet states. In this case, the BW-theory frame- 
work predicts a magnetic dipole A coefficient for the 1 P\ 
state of -32.7(2) MHz, whereas the experimental value is 
-3.4(4) MHz [5a]. Since Sg is determined mainly by the 
properties of the 3 Pi state, it is not surprising that the 
theoretical and experimental values are in good agree- 
ment. Conversely, the lifetime of the 3 Po state depends 
nearly equally on the 1 P\ and 3 Pi characteristics, so the 
lifetime prediction deserves further investigation. 

A modified BW (MBW) theory @ HE [H was at- 
tempted to incorporate the singlet data and eliminate 
such discrepancies. In this case 1 Pi, 3 Pi, and 3 P2 hfs are 
all used in the calculation, and two scaling factors are 
introduced to account for differences between singlet and 
triplet radial wavefunctions when determining the HFI 
mixing coefficients (note that 70 is not affected by this 
modification). This method has been shown to be suc- 
cessful in the case of heavier systems such as neutral Hg 
[3. We find a = 2.56(1) x 10~ 4 and f3 = -5.5(1) x 10~ 6 , 
resulting in Sgfio = -117.9(5) Hz/G and r 3p ° = 110(1) s. 
Here, the agreement with experiment is fair, but the un- 
certainties in experimental parameters used for the the- 
ory are too small to explain the discrepancy. 

Alternatively, we note that in Eq. [51 8g depends 
strongly on a^a and only weakly (< 1%) on 0ofii such 
that our measurement can be used to tightly constrain 
ao = 2.35(1) x 10~ 4 , and then use only the triplet hfs data 
to calculate /3o in the MBW theory framework. In this 
way we find /3 = -3.2(1) x 10" 6 , yielding t p ° = 182(5)s. 
The resulting 1 P 1 hfs A coefficient is -15.9(5) MHz, 
which is an improvement compared to the standard BW 
calculation. The inability of the BW and MBW theory to 
simultaneously predict the singlet and triplet properties 
seems to suggest that the theory is inadequate for 87 Sr. 
A second possibility is a measurement error of some of 
the hfs coefficients, or the ground state g-factor. The 
triplet hfs is well resolved and has been confirmed with 
high accuracy in a number of measurements. An error in 
the ground state g-factor measurement at the 10% level 
is unlikely, but it can be tested in future measurements 



TABLE II: Theoretical estimates of Sg and r p ° for 87 Sr 
Values used in Calculation 



a = 0.9996 f3 = -0.0286(3) 



Calc. 


ao 

xlO 4 


A> 

xlO 6 


(s) 


Sgfio 
m F (Hz/G) 


A 1 Pi 
(MHz) 


BW 
MBW I 
MBW II 
Ref [40J 
Ref [41, 59] 
Ref f8, 9] 


2.37(1) 
2.56(1) 
2.35(1) 

2.9(3) 


-4.12(1) 
-5.5(1) 
-3.2(1) 

-4.7(7) 


152(2) 
110(1) 
182(5) 

132 
110(30) 

159 


-109.1(1) 
-117.9(5) 
-108.4(4) b 

-130(15) c 
106 d 


-32.7(2) 
-3.4(4) a 
-15.9(5) 


a Experimental value 55] 

Experimental value from this work 
c Calculated using Eq. [6] 
d Sign inferred from Figure 1 in Ref. || 
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by calibrating the field in an independent way so that 
both gj and 5g can be measured. On the other hand, 
the 1 Pi hfs measurement has only been performed once 
using level crossing techniques, and is complicated by the 
fact that the structure is not resolved, and that the 88 Sr 
transition dominates the spectrum for naturally abun- 
dant samples. Present 87 Sr cooling experiments could be 
used to provide an improved measurement of the l P\ data 
to check whether this is the origin of the discrepancy. 

Although one can presumably predict the lifetime with 
a few percent accuracy (based on uncertainties in the 
experimental data), the large model-dependent spread 
in values introduces significant additional uncertainty. 
Based on the calculations above (and many other similar 
ones) and our experimental data, the predicted lifetime is 
145(40) s. A direct measurement of the natural lifetime 
would be ideal, as has been done in similar studies with 
trapped ion systems such as In + (3^1 and Al + [57| or neu- 
tral atoms where the lifetime is shorter, but for Sr this 
type of experiment is difficult due to trap lifetime limi- 
tations, and the measurement accuracy would be limited 
by blackbody quenching of the 3 Po state [58j . 

Table [TT] summarizes the calculations of Sg and r p ° 
discussed here including the HFI mixing parameters ao 
and Bp. Other recent calculations based on the BW the- 
ory [g, @, ab initio relativistic many body calculations 
[iq ]. and an effective core calculation are given for 
comparison, with error bars shown when available. 



VI. IMPLICATIONS FOR THE 87 SR LATTICE 
CLOCK 

In the previous sections, the magnitude of relevant 
magnetic and Stark shifts has been discussed. Briefly, we 
will discuss straightforward methods to reduce or elim- 
inate the effects of the field sensitivities. To eliminate 
linear Zeeman and vector light shifts the obvious path is 
to use resolved sublevels and average out the effects by al- 
ternating between measurements of levels with the same 
\mp |. Figure[5]shows an example of a spin-polarized mea- 
surement using the mp = ±9/2 states for cancelation of 
the Zeeman and vector shifts. To polarize the sample, 
we optically pump the atoms using a weak beam reso- 
nant with the 1 S'o- 3 Pi (F = 7/2) transition. The beam 
is co-aligned with the lattice and clock laser and linearly 
polarized along the lattice polarization axis (9 = 0), re- 
sulting in optical pumping to the stretched (mp — 9/2) 
states. Spectroscopy with (blue) and without (red) the 
polarizing step shows the efficiency of the optical pump- 
ing as the population in the stretched states is dramati- 
cally increased while excitations from other sublevels are 
not visible. Alternate schemes have been demonstrated 
elsewhere @, [26| where the population is pumped into a 
single m F = ±9/2 state using the 1 S - 3 Pi {F = 9/2) 
transition. In our system, we have found the method 
shown here to be more efficient in terms of atom number 
in the final state and state purity. The highly efficient 
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FIG. 9: (color online) The effect of optical pumping via the 
Pi (F — 7/2) state is shown via direct spectroscopy with 9 = 
0. The red data shows the spectrum without the polarizing 
light for a field of 0.27 G. With the polarizing step added 
to the spectroscopy sequence the blue spectum is observed. 
Even with the loss of ~ 15% of the total atom number due to 
the polarizing laser, the signal size of the rriF = ±9/2 states 
is increased by more than a factor of 4. 

optical pumping and high spectral resolution should al- 
low clock operation with a bias field of less than 300 mG 
for a 10 Hz feature while keeping line pulling effects due 
to the presence of the other sublevels below 10~ 17 . The 
corresponding second order Zeeman shift for such a field 
is only ~21 mHz, and hence knowledge of the magnetic 
field at the 10% level is sufficient to control the effect 
below 10" 17 . With the high accuracy ^-measurement 
reported here, real time magnetic field calibration at the 
level of a few percent is trivial. For spin-polarized sam- 
ples, a new magic wavelength can be determined for the 
mjr-pair, and the effect of the tensor shift will only be to 
modify the cancelation wavelength by at most a few pi- 
cometers if a different set of sublevels are employed. With 
spin-polarized samples, the sensitivity to both magnetic 
and optical fields (including hyperpolarizability effects) 
should not prevent the clock accuracy from reaching be- 
low 10" 17 . 

Initial concerns that nuclear spin effects would limit 
the obtainable accuracy of a lattice clock have prompted 
a number of recent proposals to use bosonic isotopes 
in combination with external field induced state mixing 
[j~7l . [H, HI H HI] to replace the mixing provided natu- 
rally by the nuclear spin. In these schemes, however, the 
simplicity of a hyperfine-free system comes at the cost 
of additional accuracy concerns as the mixing fields also 
shift the clock states. The magnitudes of the shifts de- 
pend on the species, mixing mechanism, and achievable 
spectral resolution in a given system. As an example, 
we discuss the magnetic field induced mixing scheme [20( 
which was the first to be experimentally demonstrated 
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for Yb [19] and Sr [45]. For a 10 Hz 88 Sr resonance (i.e. 
the linewidth used in this work), the required magnetic 
and optical fields (set to minimize the total frequency 
shift) result in a second order Zeeman shift of —30 Hz 
and an ac Stark shift from the probe laser of —36 Hz. 
For the same transition width, using spin-polarized 87 Sr, 
the second order Zeeman shift is less than —20 mHz for 
the situation in Fig. [9J and the ac Stark shift is less than 
1 mHz. Although the nuclear-spin-induced case requires 
a short spin-polarizing stage and averaging between two 
sublevels, this is preferable to the bosonic isotope, where 
the mixing fields must be calibrated and monitored at the 
10~ 5 level to reach below 10~ 17 . Other practical concerns 
may make the external mixing schemes favorable, if for 
example isotopes with nuclear spin are not readily avail- 
able for the species of interest. In a lattice clock with 
atom-shot noise limited performance, the stability could 
be improved, at the cost of accuracy, by switching to a 
bosonic isotope with larger natural abundance. 

In conclusion we have presented a detailed experimen- 



tal and theoretical study of the nuclear spin effects in op- 
tical lattice clocks. A perturbative approach for describ- 
ing the state mixing and magnetic sensitivity of the clock 
states was given for a general alkaline-earth(-like) system, 
with 87 Sr used as an example. Relevant Stark shifts from 
the optical lattice were also discussed. We described in 
detail our sign-sensitive measurement of the differential 
g-factor of the 1 Sq- 3 Po clock transition in 87 Sr, yield- 
ing noSg — — 108.4(4)mF Hz/G, as well as upper limit for 
the differential and exited state tensor shift coefficients 
An T = 0.02 Hz/(U T /E R ) and = 0.01 Rz/(U T /E R ). 
We have demonstrated a polarizing scheme which should 
allow control of the nuclear spin related effects in the 87 Sr 
lattice clock to well below 10~ 17 . 

We thank T. Ido for help during the early stages of 
the g-factor measurement, and G. K. Campbell and A. 
Pe'er for careful reading of the manuscript. This work 
was supported by ONR, NIST, and NSF. Andrew Lud- 
low acknowledges support from NSF-IGERT through the 
OSEP program at the University of Colorado. 
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VII. APPENDIX 

The appendix is organized as follows, in the first sec- 
tion we briefly describe calculation of the mixing coeffi- 
cients needed to estimate the effects discussed in the main 
text. We also include relevant Zeeman matrix elements. 
In the second section we describe a perturbative treat- 
ment of the magnetic field on the hyperfine-mixed 3 Po 
state, resulting in a Breit-Rabi like formula for the clock 
transition. In the final section we solve the more general 
case and treat the magnetic field and hyperfine interac- 
tion simultaneously, which is necessary to calculate the 
sensitivity of the 1 P 1 , 3 Pi and 3 P2 states. 



A. State mixing coefficients and Zeeman elements 

The intermediate coupling coefficients a and (3 are typ- 
ically calculated from measured lifetimes and transition 
frequencies of the 1 P\ and 3 Pi states and a normalization 
constraint, resulting in 



a 



'-Pi 



'Pi 



a 2 + /3 2 = 1. 



(15) 



The HFI mixing coefficients ao, (3o, and 70 are due to 
the interaction between the pure 3 Pq state and the spin- 
orbit mixed states in Eq. [I] having the same total angular 
momentum F. They are defined as 



0o = 



70 = 



( 3 P 1; P 


= i\h a \ 3 pS,f 


= 1) 




U 3p (] - 1/3 Pl 






= I\H A \ 3 P§,F 


= 1) 








{ 3 P2,F 


= I\H Q \ 3 P°,F 


= 1) 


U3p Q - U3p 2 



(16) 



Where Ha and Hq are the magnetic dipole and electric 
quadrupole contributions of the hyperfine Hamiltonian. 
A standard technique for calculating the matrix elements 
is to relate unknown radial contributions of the wavefunc- 
tions to the measured hyperfine magnetic dipole (A) and 
electric quadrupole (Q) coefficients. Calculation of the 
matrix elements using BW theory [H, [H, [H, 0, [5f| can 
be performed using the measured hyperfine splitting of 
the triplet state along with matrix elements provided in 
pMj . Inclusion of the 1 P\ data (and an accurate predic- 
tion of p ) requires a modified BW theory [H, [56| 
where the relation between the measured hyperfine split- 
ting and the radial components is more complex but man- 
ageable if the splitting data for all of the states in the 
nsnp manifold are available. A thorough discussion of 
the two theories is provided in Refs. [4^.[55|. 

Zeeman matrix elements for singlet and triplet states in 
the nsnp configuration have been calculated in Ref. [24| • 
Table IITT1 summarizes those elements relevant to the work 
here, where the results have been simplified by using the 
electronic quantum numbers for the alkaline-earth case, 



but leaving the nuclear spin quantum number general 
for simple application to different species. Note that the 
results include the application of our sign convention in 
Eq. [5] which differs from that in Ref. [24| . 



B. Magnetic field as a perturbation 

To determine the magnetic sensitivity of the 3 Pq state 
due to the hyperfine interaction with the 3 Pi and l P\ 
states, we first use a perturbative approach to add the 
Zeeman interaction as a correction to the | 3 Po) state in 
Eq. [3] The resulting matrix elements depend on spin- 
orbit and hyperfine mixing coefficients a, (3, ao, (3q, and 
7o- For the 3 Po state, diagonal elements to first order in 
ao and /3q are relevant, while for 1 P\ and 3 Pi, the contri- 
bution of the hyperfine mixing to the diagonal elements 
can be ignored. All off-diagonal terms of order /3 2 , a^a, 
ctQ0, ctg, and smaller can be neglected. Due to the selec- 
tion rules for pure (LS) states, the only contributions of 
the 3 P2 hyperfine mixing are of order ao7o, 7o\ and /3o7o- 
Thus the state can be ignored and the Zeeman interac- 
tion matrix M z between atomic P states can be described 
in the {I 1 Pi, F, mp), \ 3 Po,F,m F ), \ 3 Pi,F,m F )} basis as 



M z = 



M 



Pn 



M 1 



Pi 



\ 



V3 Pa M 3 



Pi 



V 



M. 



Po 

3 P1 



(17) 



V3 Pl J 



where we define diagonal elements as 



V3 Pq =vlpM Z Pl\HzfPl) 



+ 2(aa ~ /3/3o)< 3 Pi°,F = I\H Z \ A PH) 
U3 Pl = v 3Pl + -£(a*( 3 P?,F'\H z \ s P?,F>) 

F' 

+p 2 ( 1 P?,F'\H z \ 1 Pl\F')) 
vi Pl = v?p 1 + J2(*( 1 Pi°,F'\H z \ 1 P 1 ,F') 

F' 

+f3 2 ( 3 P?,F'\Hz\ 3 P?,F')) . 
Off diagonal elements are given by 



|3 r>0\ 



(18) 



3 p 



l=M 3 £ =a r£\{3P°,F'\H z \3P°,F)\* 



mI;i = m:;i =(3IJ2\( 3 p8,f\h z \ 3 p?,f>)\ 



(19) 
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TABLE III: Zeeman Matrix Elements for Pure ( 2 



Relevant Elements for the Pq State: 



( J Pg,F = I\Hz\ 3 P$,F = /)= -grnip^B 

< 3 P °,F = I\H Z \ 3 P?,F' = I) ={g a - gi )m FfMj B^^ 



(*P§,F = /|ff a | s P?,F' =7 + 1) =(ff. - gO^gA/ V+^iy-x? 



< 3 P °,F = I\H Z \ 3 P?,F' =1-1) ={g 3 - g.)W (f 



Relevant Diagonal Elements within 3 Pi Manifold: 



( 3 P°,F = 7|iT z | a P?,.F = /)= ( w+ *-ffff+i) +1? ~ 2? ) m^oB 
fPf.F = I + 1\H Z \ 3 P?, F = I + 1>= f 9 ' + 4i;i) 9jJ ) m F/ i B 
< 3 P°,F = 7 - 1\H Z \ 3 P?,F = J - 1>= - flWM+a ) rriFHoB 



Relevant Diagonal Elements within Pi Manifold: 



( 1 P?,F = I\H Z M,F = /)= ( a '- 9 %^ ) - 1) ) m F ^B 
( 1 P?,F = I + 1\H Z \ 1 P?,F = I + 1>= ( i f 7 ^r) m FM0 B 
( 1 P 1 °, F = 7 - l|Jf z rPi°,P = / - 1>= (- gi+ % (/+1) ) m FM0 P 



The eigenvalues of Eq. [T7] can be written analytically as 
three distinct cubic roots 



— arccos 
3 



2(yl + 3// 2 ) 3 / 2 



3 



(20) 



— arccos 



2(vl + 3^ 2 )3/ 2 



2tt 



where we have 



v =1/3 P(s + V3 Pl + vi Pl 



Ul = 



Pl\2 



v 3Pl - v 3Pl vi Pl - v; iPg vi Pl + (M 3p i) 



+(< 1 ) 2 ] 5 



(21) 



V2 = \ V3 Po V3 Pl Vl Pl 



i Pl {M 3Pl 



P M 2 



Pi 



Since the main goal is a description of the 3 Po state sen- 
sitivity, the solution can be simplified when one considers 
the relative energy spacing of the three states, and that 
elements having terms f3, a(3, and smaller are negligible 

compared to those proportional to only a. Therefore we 

in 

can ignore M ip } terms and find simplified eigenvalues 
arising only from the interaction between 3 Pi and 3 Pq 
that can be expressed as a Breit-Rabi like expression for 
the 3 Pq state given by 
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FIG. 10: (color online) Magnetic sensitivity of the 1 Pi state 
calculated with the expression in Eq.[24]using A = —3.4 MHz 
and Q = 39 MHz Note the inverted level structure. 



sensitivity of the clock transition (plotted in Fig. [5]) is de- 
termined by simply subtracting the ( 3 Pq\Hz\ 3 Pq) term 
which is common to both states. 



C. Full treatment of the HFI and magnetic field 



Vzp , mF =2 ("tPo + V3 Pi) + 2 \ V3p ° 



*Pi) 



r F . a 2 \( 3 Pg,F\H z \3P?,F>)\< 



(22) 



(V3p - V3 Pl ) 2 



For magnetic fields where the Zeeman effect is small com- 
pared to the fine-structure splitting, the result is identi- 
cal to that from Eq. [5] of the main text. The magnetic 



For a more complete treatment of the Zeeman effect 
we can relax the constraint of small fields and treat the 
hyperfine and Zeeman interactions simultaneously using 
the spin-orbit mixed states in Eq. [T]as a basis. The total 
Hamiltonian is written Ptotai = Hz + Ha + Hq including 
hyperfine Ha and quadrupole Hq effects in addition to 
the Zeeman interaction Hz defined in Eq. O of the main 
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FIG. 11: (color online) Magnetic sensitivity of the 3 Pi state 
calculated with the expression in Eq.[24]using A = —260 MHz 
and Q = -35 MHz [gg|. 



text. The Hamiltonian H to tai can be written as 



Diagonalization of the full space using Eq. [23] does not 
change the 3 Po result discussed above, even for fields 
as large as 10 4 G. This is not surprising since the 3 Po 
state has only one F level, and is therefore only af- 
fected by the hyperfine interaction through state mix- 
ing which was already accounted for in the previous cal- 
culation. Alternatively, for an accurate description of 
the x Pi, 3 Pi and 3 P 2 states, Eg. |2"31 must be used. For 
an alkaline-earth 2S+1 L\ state in the \I, J, F, TOf) basis 
we find an analytical expression for the field dependence 
of the F = I, I ± 1 states and sublevels. The solution 
is identical to Eq. [201 except we replace the frequencies 
in Eq. [21] with those in Eq. [24] We define the relative 
strengths of magnetic, hyperfine, and quadrupole inter- 
actions with respect to an effective hyperfine-quadrupole 



coupling constant Waq = A + 
X A = 



3Q 



4/(1-27) 



as X B r = 



and X Q = I(1 _ 2 y )WAQ , respectively. The so- 
lution is a generalization of the Breit-Rabi formula [6l[ 
for the 2S+1 Li state in the two electron system with nu- 
clear spin I. The frequencies are expanded in powers of 
X B r as 



-H z +AI-J 



P J(2PJ + 1)-P7(J + 1)(J+1) ( 23 ) 



2IJ(2I- 1)(2J- 1) 



no = -2Waq 



1 + '^-m F X BR 



vi = w AQ Jx2 



. . 2 (Sefi - 9i)Xa + 3g c f{XQ 
1 H — m F X B R - 



A cff 



x v h 

cii 



X 



2 

BR 



V2 = W A q tfHj + l)X2 



X l (7$^ + Bl) + X% 3(1 - 2/ 1 ) 6 (3+2/) (j^y - 9I ) - X A X Q ( 9cff (2 - + 9l 



(24) 



x v l 

Ctt 



m F X BR 



x * 7(7+T) + X Q 3 ^ " /(/ + i)(g cff +g^ ) vS 



f ((gcff + gj) 2 -(g/"lF) 2 ) 



X BR 



1(1+1) 



x v l 

err 



-m F Xf 



with abbreviations 



= T V + i)(x A + ^- i(i + i)x Q (x A - 1)) - 1 



X e £ =X e s [ XqX c b + ( X A - X, 



r2 v2 3(3 + 2/)(l-27) 
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X e g =Xa + X ( 



(3 + 2/)(l-2J) 



g eS = <*i±?A + (i(L + 1} _ s(5 + 1)} _ 



(25) 



The resulting Zeeman splitting of the 5s5p 1 P\ and 
5s5p 3 Pi hyperfine states in 87 Sr is shown in Fig. [TOl and 
Fig. [IT] For the more complex structure of 3 P2, we 
have solved Eq. [23] numerically, with the results shown in 
Fig. HH The solution for the 1 Pi state depends strongly 
on the quadrupole (Q) term in the Hamiltonian, while 
for the 3 P\ and 3 P2 states the magnetic dipole (A) term 
is dominant. 
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FIG. 12: (color online) Magnetic sensitivity of the 3 Pi state 
calculated numerically with Eq. (23 using A=-212 MHz and 
Q=67 MHz [12. 



